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Abstract: 

We find new representations for Itzykson-Zuber like angular integrals for arbitrary 
/?, in particular for the orthogonal group 0{n), the unitary group U{n) and the sym- 
plectic group Sp{2n). We rewrite the Haar measure integral, as a fiat Lebesge measure 
integral, and we deduce some recursion formula on n. The same methods gives also 
the Shatashvili's type moments. Finally we prove that, in agreement with Brezin and 
Hikami's observation, the angular integrals are linear combinations of exponentials 
whose coefficients are polynomials in the reduced variables (xj — Xj){yi — Uj). 



1 Introduction 

What we call angular integral [2S] is an integral over a compact Lie group Gp^n- 

Gi,2,n = 0{n) , G'i,„ = f/(n) , G2,n = Sp{2n) (1-1) 

of the form: 

W{X,Y)= [ dOe^^°^°-' (1-2) 

where X and Y are two given matrices, and dO is the Haar invariant measure on the 
group. We shall also extend Ip^n to arbitrary [3 (Notice that our [3 is half the one most 
commonly used in matrix models, for instance we have /? = 1 in the unitary case). 

In this paper we are going to consider the case where X and Y are diagonal matrices, 
however, let us first recall the Harish- Chandra case. 
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Harish- Chandra case 

In the case where X and Y are in the Lie algebra of the group [22] (i.e. real anti- 
symmetric in the 0{n) case, anti-hermitian in the U{n) case, and quaternion- ant i- self- 
dual in the Sp{2n) case), the angular integral can be computed with Weyl-character 
formula, and is given by the famous Harish- Chandra formula [10] (which is also a 
special case of the Duistermaat-Heckman localization [8]): 



where C is a normalization constant, w runs over elements of the Weyl group, and 
the generalized Vandermonde determinant A^(X) is the product of scalar products of 
positive roots with X (see [IDl EHl ES] for details). 

Diagonal case 

However, for applications to many physics problems [S1I2S], it would be more interesting 
to have X and Y in other representations, and in particular X and Y diagonal 
matrices. 

Since a antihermitian matrix is, up to a multiplication by i, a hermitian matrix, 
and since every hermitian matrix can be diagonalized with a unitary conjugation, for 
the unitary group, the Harish- Chandra formula applies as well to the case where X 
and Y are diagonal, this is known as Itzykson-Zuber formula |18j : 



where A(X) = Ai(X) = H ( the usual Vandermonde determinant. 

For the other groups, computing angular integrals has remained an important chal- 
lenge in mathematical physics for a rather long time. Many progresses and formulae 
have been found, however, a formula as compact and convenient as Harish-Chandra 
is still missing. And in particular a formula which would allow to compute multiple 
matrix integrals, generalizing the method of Mehta [26] is still missing. 

Calogero Hamiltonian 

It is known that, in the diagonal case, satisfies the Calogero-Moser equation [3], 
i.e. is an eigenfunction of the Calogero hamiltonian: 



(X, Y) G Lie algebra 




(1-3) 






(1-5) 
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Many approaches towards computing angular integrals have used that differential 
equation. A basis of eigenfunctions of the Calogero hamiltonian is the Hi-Jack poly- 
nomials [ll[2l[71[2l]. 

In particular remarkable progress in the computation of was achieved recently 
by Brezin and Hikami [3]. By decomposing „ on the suitable basis of Zonal polyno- 
mials, they were able to find a recursive algorithm to compute the terms in some power 
series expansion of I^^n, and they obtained a remarkable structure. In particular they 
observed that the power series reduces to a polynomial when P E N. 

Morozov and Shatashvili's formulae 

Another important question for physical applications, is not only to compute the an- 
gular integral (the partition function in statistical physics language), but also all its 
moments, for instance: 

M,,,= [ dOe^^^^^^-^ ||0,,|r (1-7) 

and more generally for any indices ^l, . . . , i2k,ji, ■ ■ ■ ,j2k- 

TrXOYO-^ 0^^^n.^^...Oi,, ,i,,Oi\ 0;\ ...O-^ , (1-8) 

In the U{n) case f3 = 1, Morozov [28] found a beautiful formula for Mj j, and 
Shatashvili [30] found a more general formula for any moments of type 11-81 using the 
action-angle variables of Gelfand-Tseytlin corresponding to the integrable structure of 
this integral. 

For /? = 1/2,1,2, in the Harish- Chandra case where X and Y are in the Lie algebra, 
a formula for all possible moments was also derived in p^, generalizing Morozov's 

mm- 

In this article we shall propose new formulae for Mij in the diagonal case for 
arbitrary f3, and our method can also be generalized to all moments. 

Outline of the article 

• Section 1 is an introduction, and we present a summary of the main results of 
this article. 

• In section ?? we setup the notations, and we review some known examples. 
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In section H] we show how to transform the angular integral with a Haar measure 
into a flat Lebesgue measure integral on a hyperplane. From it, we deduce a 
recursion formula, as well as a duality formula (the angular integral is an eigen- 
function of kernel which is the Cauchy determinant to the power j3). 

In section [5], we discuss the moments of the angular integral. We show that 
moments can be obtained also with Lebesgue measure integrals, and we show 
that they satisfy linear Dunkl-like equations. This can be used as a way to 
recover Calogero equation for the angular integral. 

In section [HI we rewrite the angular integral as a symmetric sum of exponentials 
with polynomial prefactors. Those polynomials are called principal terms, and 
can be computed recursively. In particular, we prove the conjecture of Brezin and 
Hikami [3] that the principal terms are polynomials in some reduced variables 
{xi- Xj){yi-yj). 



In section 16. 2[ we prove a formula for n = 3 in terms of Bessel polynomials, and 
we propose a conjecture formula for arbitrary (3 and arbitrary n. 



In section 16.41 we focus on the symplectic case P = 2, for which we can improve 
the recursion formula. 

Section [7] is the conclusion. 

Appendices contain useful lemmas, and proofs of the most technical theorems. 



.1 Summary of the main results presented in this article 

• We rewrite the angular integral with the Haar measure on the Lie group Gjs^n, 
as a flat Lebesgue measure integral on its Lie algebra (notations are explained in 
section Hj): 



,TvS 



dS 



m=idetiS-y,Xy 



1-9) 



as well as its moments: 



Mij = J dO \\0^,,\\^ e 



2 ^TrXOYO-^ 



:i-io) 



We show that the Mjj's satisfy a hnear functional equation (very similar to Dunkl 
operators): 

^hj, -fT^ + f^l^ — _ — - = Mijyj (1-11) 

which imphes the Calogero equation for /^^^i = X]i = X]j 

E^+''E^(l^-§^) = £^?)/.,. (1-12) 

Moreover, the integral of eq. (Il-iop . is a solution of the linear functional equation 
eq. fll-lip for any choice of integration domain (as long as there is no boundary 
term when one integrates by parts). We thus have a large set of solutions of the 
linear equation, and also of Calogero equation. 

We deduce a duality formula: 

/,,.(X; Y) = det{Xy-^ [ dX, . . . rfA„ AiA)'^ u^n^n^ ^ (1-13) 

and a recursion formula: 



^An(X;F) 



;i-i4) 



similar to that of [IHl [16] . 



For /3 G N, the solution of the recursion can be written in terms of principal 
terms: 

En 

Il3,n{X, Y)=2_^ A(X)2/3A(r„)2/3 ^'3'"(^' (^"^^^ 

where is the sum over all permutations. 

The recursion relation eq. (11-141) can be rewritten as a recursion for the principal 
terms Xg^n(X, Y): 



a.i=yi 



(1-16) 
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For general n and (3 integer, we prove the conjecture of Brezin and Hikami [3], 
that the principal term Xg Y) is a symmetric polynomial of degree (3 in the 



Tij variables, 



(x, - Xj){yi - Vj) 
2 



r,,, = - ^ ' (1-17) 



In the case n = 3 we find this polynomial explicitly for any (3 (for (3 integer the 
sum is finite): 



QXiyi+X2y2+xsy3 Y(d — k) -i-r (h') 1 



where ym is the m*^ Bessel polynomial, i.e. the modified Bessel function of the 
second kind (see definition of 3^/3_i below in eq. 03-41) ). 

In the case (3 = 2 (i.e. symplectic group Sp{2n)), the recursion relation for the 
principal term can be written: 



^2,n = ]X^=i{Xi - Xn) iVi - Vn)^ 



ai=yi 



A(Y \^/\(V\^ det +B+dY] /,^„_,(X„_i;y„_i) 

^{^nj ^{J^nJ det(X„-x„) A(X„_i)2 A(r„_i)2 



:i-i9) 



and -B is the antisymmetric matrix Bij = , = 0, and dy = 

diag(9y^, . . . , <9y^_J. In section 16.41 we propose an operator formalism to com- 
pute it, and we propose a conjecture formula in terms of decomposition into 
triangles. 

2 Definitions and examples 

3 secdefex 

3.1 Notations for angular integrals 

Let X and Y be two diagonal matrices of size n: 

X = diag(xi, . . . , x„) , Y = diag(?/i, ...,?/„) (3-1) 
We define the angular integral: 

/ft„(xi,...,x„;yi,...,i/„) = f c^Oe^^o^^"' (3-2) 
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where Gp^n denotes one of the Lie groups: 

Gi/2,n = 0{n) , Gi,n = U{n) , G2,n = Sp{2n) (3-3) 

and dO is the invariant Haar measure on the corresponding compact Lie group. 
We will later extend those notions to arbitrary values of (3. 



3.2 Bessel polynomials 

For further use, we need to introduce some Bessel functions [H [231 ED E21 E]- Those 
special functions are going to play a major role throughough this article. 

The Bessel polynomials (see [221 [32]) 3^m(a;) are defined by: 

y-^i-) - ± ^^:'^!{^ i-i^f - v^) (3-4) 



where /C is the modified Bessel function of the second kind [H [31]. J^m is a polynomial 
of degree m when m is an integer: 

= 1 , y^ = x+l, 3^2 = + 3x + 1 , ^5^3 ^ ^5^2 + gx + 1 , etc . . . (3-5) 

They satisfy: 

x'^yl + (2x + 2)X^ - m{m + 1)3^^ = (3-6) 

We shall also need: 

which is a polynomial of degree /3 — j if /? is an integer. 

In particular Qp^ is the Carlitz polynomial O [32] and is closely related to yp-i- 

Qpfi{x) = x^yp-ii-) = \^ e^ x^^-^ K^_.{x) (3-8) 
X V vr ^2 



satisfying: 
The first fews are: 



x^Q'' - 2x(/5 + x)g' + 2/?(x + l)g^,o = (3-9) 



Qi,o = 2; ! ^2,0 = 3^^ + 3;, Qs^o = 2:^^ + 3x2 + 3x , (^4,0 = a;^ + 6x^ + 15x2 + 15x , etc... 

(3-10) 
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For higher j, the Qp/s are derivatives of Bessel polynomials: 
They satisfy: 



g^,,(i/x) = 2^- x^-f' — y^_^{x) = T x^-p y'^^l.ix) (3-11) 



- X Qp,^ = \ Q^,,+i + iQ^,,- + (j - + /3 - l)g^_^._i (3-12) 

Qp^^_,,^2{P-j-x^)Qp,j (3-13) 

The first fews are: 

Q2,i = 2x , Q3,i = + 12x , ^4,1 = 12x^ + 60x^ + 90x , (3-14) 

Q3,2 = 24x , Q4,2 = 120x2 + 360x , ^4,3 = 720x , etc . . . (3-15) 

3.3 Examples angular integrals with n = 1,2,3 

• n — 1: The n — 1 case needs no computation, and gives: 

If,,i{x;y) = e^y (3-16) 

• n — 2: The n — 2 case requires a little bit of easy computation, and it has been 
known for some time, we have (this formula is rederived in this article): 

ip,2{x,Y) = y^-i(i/T)+ ^_^^p y/3-i(-i/r) 

QXiyi+X2y2 Qxiy2+X2yi 

= — — <9/3,o(r) + <3/3,o(-r) (3-17) 

where ^ 

^ = "2(^1 -^2)(l/i -1/2) (3-18) 

It can also be written in terms of the modified Bessel function I: 

„ i (a; i+a;2) (2/1+2/2) 

b,2ix;y) = — — v^(^) ' (3-19) 



where 



^W = (V2)-E „r(!: + r^l) . ^™=^: + i^Il (3-20) 



n — 3: 
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we show in this article that (proof in appendix [C|) : 



(3-21) 

where 

r,, = - - ^^^^y^ - (3-22) 
and +perm. means that we have to symmetrize over all permutations of the |/j's. 

• ri > 3: We show in this article that for arbitrary n and /3, the angular integral is of 
the form conjectured by Brezin and Hikami: 

(A(x)A(y))2/3 ^^'"^'"^^'^ +perm. (3-23) 
where „( polynomial in the Tj j's, and for which we write a recursion relation. 

4 Transformation of the angular integral 

In this section, we transform the Haar measure group integral into a flat Lebesgue 
measure integral. 

4.1 Lagrange multipliers 

For (3 = 1/2, 1, 2, an element O G Gp^n is an orthonormal basis, i.e. a collection of n 
orthonormal vectors Ci, . . . , e^, whose coordinates Oij = (ej)j are of the form: 

2/3-1 

(e.), = = 5^ (e.)J 6„ (4-1) 

where the e^'s form a basis of a Clifford algebra (indeed this reproduces the three 
groups Gf3,n for P = 1/2, 1,2): 

eo = 1 , ej = eo , Va > : = -1 , e)j, = -e„ , Ca-Ca' = -e^'-ea 

(4-2) 

with structure constants (only for (3 = 2): 

= '7a,a',a" (4-3) 

a" 

and where r]a,a',a" has the property useful for our purpose, that for every pair {a, a'), 
there is exactly only one a" such that ria,a',a" 7^ 0. In particular ria^a,a" = Sa",o- 
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The basis must be orthonormal, i.e. 



20-1 



fc=l k=l a,a'=0 

We introduce Lagrange multipliers to enforce those orthonormality relations 

5{ei.e\ - 1) = y dSi,i eS^^^^^-^^,c.i{e.)t? (4.5) 

and if i < j: 

= y dSij e-2^fc^^'^(("')^)(("^-)fc) (4-6) 

where each integral is over the imaginary axis. 

Since the scalar product is invariant under group transformations (i.e. change of 
orthogonal basis), the following measure is invariant and thus must be proportional to 
the Haar measure: 

docxji d(e,)° n ^^(^-4 - 1) n '^(^-^1) ('^-^) 

i,j,a i i<j 

i.e. 

where 

n 2/3-1 

= n ^"^^.^ n = n ^-^^ n n ^-^^^ (^-9) 

i i<j i=l i<j a=0 

is the Gi3^n invariant measure on the space Ep^n' 

{Ei/2,n — {n X n real symmetric matrices} 
Eiji = {n X n hermitian matrices} (4-10) 
E2,n — {n X n quaternion self — dual matrices} 

where we have completed 5" by self duahty (5" = S^) : 

Sl^Sl^ ,andVa>0 S^.^-S^j (4-11) 
Therefore we have (up to a multiplicative constant): 

WX;y) OC /.s/de....de„ei:.»«ei:..-K«.).P 

10 



(4 - 12) 

The integral over the (ei)^'s is now gaussian and can be performed. The gaussian 
integrals for each k are independent. 

The quadratic form in the exponential is, for each k: 

Yl 'Yi(^i)ki(^j)k'Va,a',a"iSijXiykSa,>>^o- S^'-) (4-13) 
a,a',a" i,j 

If we define the vector Vk = (fi,fc, . . . , Vn,k) where Vi^k = Sa(ei)fceL we have to compute 
the gaussian integral: 

dvk e-^fe(^-^^^)^'= (4-14) 



For the 3 values of f3 = 1/2, 1, 2, this integral is worth: 

where det is the product of singular values (see [25]). 
Thus we get the following theorem: 

Theorem 4.1 The angular integral /^^„(X;y) is also equal to the following flat 
Lebesgue measure integral: 



Tr 5 r g Tr 



/ Ir D r 
llfc=idet(^-?/fcX)^ J llfc=idet(^-?/fc 



rf5 



(4-16) 

In the last formula we have made the change of variable S X^^'^SX^^'^. Also, the 
integration domain for S, which was iEjj^n before exchanging the integrations over S 
and e, is now shifted to the right, so that all singular values of {S — ykX) have positive 
real part. The integration domain for S can be deformed such that the integral remains 
convergent and the integration path goes to the right of all zeroes of the denominator. If 
(3 is half-integer or integer, the denominator is not singular near oo, and the integration 
contour can be closed. This will be made more precise below. 

Remark 1: For the moment, this formula holds only for f3 = 1/2, 1,2. Later we 
will extend it to other values of /3. 

Remark 2: Another remark, is that a similar formula can be obtained by exchang- 
ing the roles of X and Y. 
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4.2 Duality formula 

Notice that the matrix 5* itself can be diagonahzed with a G/3,n conjugation: 

S = 0A0-' , A = diag(Ai, . . . , A„) , O G 
and the measure dS is up to a constant [25] : 

dS oc dO dA A (A) 2^ 
Therefore, the angular integral reappears in the RHS: 

,Tr SX 



(4-17) 



(4-18) 



V(X;y) cx det{Xy-^ [ dS 

J 11a:=1 



oc det(X)i"^ / f/A A(A)2/5 .^-^^'"^^ 

y nLin-=i 



(A, - Z/fc)^ 



(4 - 19) 



Here, if we assume that Vi, Xi G M"*", the integration contours for the Aj's are of the 
form r + where r > max(Re?/fe). If /5 is integer or half integer, the denominator in 
the integrand is not singular near oo, and the integration contour can be closed. Thus, 
if 2(3 is an integer, the integration contours for the Aj's can be choosen as circles of 
radius > max(|?/fc|). 

This equation looks better if we rewrite it in term of the Cauchy determinant 

A(X)A(F) 



D„(X,y) =det 



1 



and the rescaled function 

Ip.n{xi, . . . ,Xn;yi, . . .,yn) 
We then have: 



(A(X)A(r))^ 



(4-20) 



(4-21) 



Theorem 4.2 the rescaled function Ip^n satisfies the duality formula: 




(4-22) 



i.e. Ii3^n{X]Y) is an eigenfunction of the kernel D!^. 



Remark 1: The duality formula above was derived for f3 = 1/2,1,2, but it makes 
sense for any j3. 
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Remark 2: It is easy to check that this relation is satisfied for the Itzykson-Zuber 
case (3=1, indeed in that case we have Ji^„(X;y) = det(e''*^^) = '}2p{~^)^ ]^. e^^^pt'), 
and: 



J dA/i,„(X,A) D^{A,Y) 

a,p 1=1 



,S/o-(i) ^p(i) 



<j,p i=l 

= n\ det(e''*^^) 

(X /i,„(X,F) (4-23) 
4.3 Recursion formula 

First, let us notice that we can always assume that Xn = 0, otherwise we perform a 
shift X X — Xn'- 



[ gTr (x-x„)oyo-i gXnTr oyo-1 

3X„Try f ^QgTV(X-x„)OyO-i ^4,24) 



Thus we define: 

X„_i = diag(a;i, . . . ,x„_i) , X = X„_i - x„ Idn_i (4-25) 



Then, we notice that the orthonormality of the basis e^: 

Ci.e] = 6ij (4-26) 

implies that if we already know ei,...,e„_i, then e„ is completely fixed (up to an 
irrelevant phase). In other words, it is sufficient to enforce only the orthonormality of 
ei, . . . , e„_i with Lagrange multipliers, i.e. introduce a matrix S of size n — 1. 

Also, because of our shift X ^ X — Xn, we notice that e„ does not appear in the 
integrand. 

Then, we write as in eq J4-81 

n— 1 n— 1 n— 1 

dO OC Yl da J] 5{ei.ei - 1) J] ^(^-^5) 

i=l 1=1 i<j=l 
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n-1 n 2/3-1 
1=1 j=l a=0 



fj^^ ^2i^i,i g~ Si Sfe g~2 ^j^^- 5jj (ei)fe (ej)|. 



/3,n-l 



oc 



n—1 „ 



(4-27) 



which imphes, after performing the gaussian integral over the ei, . . . , e„_i: 



Ip,n{X;Y) OC e^ 



.try /" 



dS 



nLidet(5-|/.X) 



oc 



/3 



Ci5 



nLidet(5-y,)^ 



(4 - 28) 
Again, 5* can be diagonahzed: 

S^OW-^ , A = diag(Ai,...,A,_i) , OeG^.n-i (4-29) 
i.e. the rank n angular integral Ip^n is expressed in terms of the rank n — 1: 



If3,n{X;Y) oc 



^Xn tr Y 



oc 



nr=i(^.-^n)^^-^ 



dA 
dA 



/^,._i(X,A)A(A)^^ 

nLinr=7(A.-i/.)^ 

/^,^_i(X„_i,A)A(A)^/^ e-""^^^^ 

nLinr=i'(A. -?/.)'' 



(4 - 30) 

Which is our main recursion formula: 

Theorem 4.3 The angular integrals //3,„(X;F) satisfy the recursion: 



OC 



ni=l i^i ^n) 



2/3- 



- J dAi, . . . o!A„_i ro^inr=-/(A,-2/. 



(4-31) 



Here again, the integration contours for the Aj's are such that the integral is convergent, 
and such that they surround all the ykS. For instance, if 2f3 is an integer, and if 
\/iXi G M^, the integration contour for the Aj's can be choosen as circles of radius 
> niax(||/fc|). 

RemcLrk 4.1 Now this recursion formula can be used to define Ijs^n for arbitrary /?, so that 
it coincides with the angular integral for /3 = 1/2, 1, 2. 
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We have also the iterated form: 

T S^" II ^n— 1 i 



n:!:a'n;tin>=i(A..'-A.+i,)/' 

(4 - 32) 

where we have defined Xnj = yj, and where the integration contours are circles such 
that: 

\Xij\ = Pi , pi> p2> ■■■> Pn-1 > max|?/j| (4-33) 



Remark 4.2 A similar recursion relation was also found Kohler and Guhr |15l [T6t 117] . 

but the authors found real integrals instead of contour integrals. The advantage of our 
formulation, is that we can easily move integration contours and find new relations, as we 
will see below. 



5 Moments of angular integrals and Calogero 

In this section, we compute moments of the angular integral, and we show that our 
formula indeed satisfies Calogero equation. 

5.1 Generalized Morozov's formula 

Define the quadratic moments (see ^28j for (3 = 1): 

Mi,= [ dO ||0,,||2e^^°^°~^ (5-1) 



gTr5 



The same calculation as above yields: 



= P Ids {{S - y.xrU 



(5-2) 



m=idet{S-y,Xy 



As a consistency check, and as a warmup exercise, let us show that this formula 
satisfies: 

/^,„ = 5^M„- (5-3) 



3 
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which comes from VO G G/j.n, J2j ll^ijlP = 1- 
We have: 



3 



det(X)i-'^ / dS / 
^An (5-4) 



Tr5X 

det(5-l/fc)'3 



Notice that this equahty holds independently of the integration domain of S, provided 
that one can integrate by parts without picking boundary terms. 

Remark: Of course a similar equation can be found by exchanging the roles of X 
and y, and one gets symmetrically: 



^M„ = /;5,„ (5-5) 



5.2 Other moments 



Since, after introducing the Lagrange multipliers, the integral becomes gaussian in the 
Oj j's, any polynomial moment can be computed using Wick's theorem. It is sufficient 
to compute the propagator: 

< Oi^uOli >=f36k,i {{S - VkX)-')^^^ (5-6) 

Then, the expectation value of any polynomial moment is obtained as the sum over all 
pairings of the product of propagators. 
For instance: 



/ dO eT>-^oyo-i 0^ O^- 

/ '-^n,n ^^12,32 ^■lz,3-J, ^t4.,34 

TrS 



dS 



5n,3^ h2,. {{S - y,, X)-i)^^_^^^ {{S - y,, X)~^) 

+^...4 ^32,3. ((5 - V3. Xy'),^ {{S - y,, X)-i) 
(5-7) 



In principle, one could compute with this method the generalization of all Shatashvili's 
moments 1301. 
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5.3 Linear equations 



In this section we prove that the Mjj's satisfy the following linear functional relations, 
which are very similar to Dunkl equations []: 



% yi- vj 



(5- 



We are going to give 2 different proofs of eq. (15-81) . The first one below is based on 
integration by parts. It can be done for the 3 groups f3 = 1/2, 1, 2, however it is rather 
tedious for (3 = 1/2 and (3 = 2, and we present the proof only for (3 = 1. Another proof 
valid for all 3 values of (3 is presented in section 15.61 below. 

Let us check that eq. (l5-2p satisfies eq. fl5-8p (for P = 1). We first rewrite: 



yi-Vj \S-yi S-yj 



J2((S-yir\m {{S-y,)-') 



(5-9) 



For (3 = 1, we may consider all variables Si^m to be independent variables, and we 
integrate by parts: 



yi - Vj 



EE^det(X)-/^ 



det(5-2/fc)/5 



-V - det(X)i-'3 I dS- ,^ 



_ 9 1 



d 



,Tr SX 



iiS-y, 



}m,i 



-/5E^det(X)^-«|d5^ 



,Tr SX 



ii^ yj) )m,iXi^i,r 



detiS-yk)" 



{{S-y,)-\m{{S-y,r')r 



,TrSX 



det{S - yk) 



i{S-yj) ^)m,m{{S-yj) ^)i,i 



dS 



det(X)^-'3 , ^ , 

' n.det(^-i/fe)/5 

f ^TrSX 

-— det(X)i-'^ / dS 



UkdetiS-VkY 



((^ 



as - y,: 
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dS 



) 



Tr {S - y,r\{S - y,r% 



(5-10) 



QED. The same computation can be repeated for /3 = 1/2 and (3 = 2, with additional 
steps because the variables Si^m are no longer independent, and also because for /3 = 2, 
det(S' — yj) is defined as the product of singular values. Another proof is given in 
section I5.6[ 

Remark: Of course a similar equation can be found by exchanging the roles of X 
and Y, and one gets the symmetric linear equation: 



Remark: again, this proves that eq. f l5-2p is solution of the differential equation 
eq. (15-81) . for any choice of integration domain provided that we can integrate by parts. 
In fact, by taking linear combinations of all possible integration contours, we get the 
general solution of the linear equation eq. (15-81) . However, a general solution of eq. (15-81) 
is not necessarily symmetric in X and Y, and does not necessarily obey eq. (15- lip . 

5.4 Calogero equation 

Here, we prove that Ip^n satisfies the Calogero equation. 
Start from the linear equation: 




(5-11) 




(5-12) 



then sum over j, using eq. (l5-3p : 




(5-13) 



Then apply -tt-: 
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dl_ dl_ 



(5 - 14) 

If we take the sum over i, using eq. (15-51) . we get: 



dl _ dl 



i.e. we recover the Calogero equation: 



H CalogcTO- 1 13, n — C^^v'j) ^I3,n 
3 



5.5 Matrix form of the linear equations 

The hnear equations, are rt^ hnear equations of order 1, for rt^ unknown funct 
They can be summarized into a matricial equation: 



MY = KM 



where i^' is a matricial operator 
d 



and more generally this implies: 



ik 



i ^ k 



MYP = K^M 



and therefore, for any polynomial P: 

M.P{Y) = P{K).M 
In particular if we choose the characteristic polynomial of Y: 

n 

= l[{y,-K) .M 

i=l 

Let us introduce the vector 



e= (1,1,...,1) 
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It is such that M is a stochastic matrix, i.e.: 



M.e = J^,„ e , e*.M = e* 
We thus have, for any polynomial P: 

e*P(K)e.J^,„ = J^,„ TrP(r) 



(5-23) 



(5-24) 



Notice that the Calogero equation is the case P{K) = K^. 

If P is the characteristic polynomial of Y we get another differential equation for 



(5-25) 




And if P{K) = n,^,(l/i - K), we get: 




(5-26) 

This last relation allows to reconstruct j if we know 

Finally, before leaving this section, we just mention that those operators Ki^j are 
also related to the Laplacian over the set of matrices Ep^n) as was noted recently by 
Zuber 



5.6 Linear equation from loop equations 

There is another way of deriving those Dunkl-like linear equations for the angular 
integrals, using loop equations of an associated 2-matrix model. 

Consider the following 2-matrix integral, where Mi and M2 are both in the E^^^n 
ensemble: 

Z= ZrfMirfMa e-T'-^^i^^i^+^^^^^^-^i^^) (5-27) 



After diagonalization of Mi = OiXO^ ^ and M2 = O2VO2 ^, we have: 

Z = j dXdYdOM A{Xf^A{Yf^ g-TV(y,(x)+y,(y)) ^Tr xo-^o,yo-^o, 

We redefine O2 = Oi.O, and the integral over Oi gives 1, and the integral over O gives 
the angular integral: 



z = J dXdY a(x)2/5a(f)2/5 ^-Trmx)+v,(Y)) /^^^(x,y) 



(5-29) 
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We can do a similar change of variable for moments: 

< Tr ^— > 

X- Ml y- M2 



- f dMidM2 e-Tr(Vi(M,)+V.(M,)-M,M,) ( _J: 1 

Z J \x- Ml y- 

i j dXdYdOM AiXf'AiYf^ g-Tr(y,(x)+y,(y)) 

gTi.xo-o.yo-Oi Tr ( Oi'02 ^—O.'Oi) 

\x-X y-Y J 

I J dXdY A{X)'^A{Y)'^ e~Trmx)+v,iY)) 

yMiJX,Y) 

' x-Xi y-yj 



(5 - 30) 



where Mij{X, Y) is the Morozov moment defined in eq J5-l[ 

Loop equations amount to say that an integral is invariant under a change of vari- 
ables. Thus, we change Mi Mi + e ^J-^^^ y-Ah "^"^(^^) order 1 in e we get 
(the loop equations for (3 = 1/2, 1, 2 ensembles can be found in several references [], the 
Jacobian is easily computed in eigenvalue representation, see appendix [Bl eq. (12-81) ): 



X- Ml y- M2/ \ X- Ml y - M2 
+(3 ( Tr Tr ^ ^ 



X- Ml X- Ml y - M2 



dx \ X- Ml y- M2 

(5-31) 

I.e., gomg to eigenvalues Mi = OiXOi^ and M2 = OzFOg"^: 

= ^ /(% -F/(x,))M,,(x,r) 

id 

+i3y^y^ \ 1 — vi — vi — ^ 

\ [x - xi)[x - x,{){y - y.j) 



{x-Xi){y-yj) 

M,,,(x,y) 

{x - Xi)'^{y - yj) 



(5 - 32) 

The last term can be integrated by parts: 
^\{x-Xif{y-yj) 
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J2 fdXdY A{Xy^A{Yy^ e 



{Vi{X)+V2(Y)) 



dxi {x - Xi) {y - Uj) 
j dXdY _ 

''iJ 

y [dXdY AiXff'AiYf^ 
{x-Xi){y-yj) 



^-'Tr{Vi(X)+V2(Y)) 9Mij{X,Y) 

dXi 



^ I Vl{x,)M,,,{X, Y) \ o . / M,,(X,F) 

f-\ {x-x,){y-yj) I ^j^^\{xi- 
j^/ 1 dMij{X,Y)\ 



,(a;-a;i)(|/-|/j) dxi / 
(5 - 33) 

Therefore we have: 

1 aMij(X,Y)\ 



E 



J ,.x-xi){x-Xi){y-yj) 



y,M,,{X, Y)_\_^^y^/ M,,(X, Y) 

M^,,{X,Y) 



^ ^ \ - ^0 (a^i - ^i) {y - Vj) 

Xi){xi - xi){y - yj) 



E 



{x - Xi){y - yj) 

1 M,,iX,Y) - M,,iX,Y) \ 

i^i J ,{x-Xi){y-yj) (xi-xi) / 



(5 - 34) 
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Since this equation must hold for any Vi and V2, x, y, i.e. for any measure on X and 
y, it must hold term by term i.e. we recover the linear equation: 



'-Mlj^ . y, M.,(A-, y)-pY, -^'^'^f 

dxi ^ {xi - xi) 



Of course, the loop equation coming from the change of variable M2 



1 1 



+ O(e^) gives the symmetric linear equation: 



9yj ' j^. ivj-yi) 



QED. 



(5-35) 
M2 + 

(5-36) 



6 Principal terms and the Tij variables 

As we mentioned in the introduction, it was noticed in particular by Brezin and Hikami 
[3], that the angular integral can be written as combinations of exponential terms, 
and polynomials (for f3 integer, series otherwise), of some reduced variables Tij = 
— I {xi — Xj){yi — Uj). Here, we show how our recursion gives such a form. 
We thus define: 



Definition 6.1 We define the principal term Ip^n{X; Y) from the recursion: 



and 



n—l 



Hynf^W{x^-Xn) Res 



dXn-l 



i=l 



\^y^ (Ai - yi)f^ ' ' ' (A„_i - Vn-iY 



(6-1) 



It is such that (the sum over permutations comes from the sum of residues at all poles 
in recursion eq. \4-31 ) of theorem. \4-3^ : 



//3,n(X,F) = 2^ A(X)2/5A(F^)2/5 l(i,n{X,Y„) 



(6-2) 



In |3], Brezin and Hikami observed and conjectured that when (3 is an integer, 
Xp^n{X, Y) is a polynomial in the variables Tij = — |(xi — Xj){yi — Uj)- 
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For instance, if /3 = 1, we have for arbitrary n: 

ii,n = Yl ^iJ (6-3) 

i<j 

And, if n = 2, we have for arbitrary /3: 

dX g(a^i-2.'2)A 

I,,2 = (x, - X,) {y, - Res ^ _ _ 

= {xi - X2) {yi - y2f'^ 



13-1 



A=0 



= {xi - X2) {yi - y2Y^ 

k\{P-l-k)\ ^""^ ""'^ (/?-!)! 

= 2^ Q/3,o(ri,2) (6-4) 

i.e. we recover the well known result that Xg^2 is the Bessel polynomial of degree /3. 

We are going to prove the conjecture of Brezin Hikami for all n and for all /5 G N, 
but first, let us prove some preliminary properties: 

Lemma 6.1 Jp^n is a polynomial in all variables Xi and yj, and it is symmetric under 
the exchange X ^ Y, and under the permutation of pairs {xi,yi) ^ {xj,yj), and under 
translations X ^ X + cte.Id or Y ^ Y + cte.Id. 

proof: 

If P is an integer, the recursion relation eq. (16-11) leads to (we write Xij = Xi — Xj, 

Vij = yi-yj)- 



i=l 



(6-5) 



l<j<n+lj^ 



A=0 



which shows by recursion, that is a rational function of all Xj's and t/j's. 

We know, from its very definition, that the angular integral 

Ii3,n{X, Y) = 2_^ _______ „(X, Y„) (6-6) 
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is symmetric in all Xj's and y^'s, and in the exchange X F. Since the exponentials are 
linearly independent on the ring of rational functions, each term must be symmetric in 
permutations of pairs (xj, ?/i)'s, i.e. 2'^^„(a;i, . . . , x„; yi, . . . , Un) is a symmetric function 
of the pairs {xi,yi)^s, and also symmetric under X ^Y. 

Moreover, X^,n+i is clearly a polynomial in the variables Xn+i and yn+i, and because 
of the symmetry, it must also be a polynomial in all variables. Translation invariance 
is also clear from the recursion formula. 

□ 

Theorem 6.1 (Conjecture of Brezin-Hikami): 

il3,n is a symmetric polynomial of degree (3 in the Ti j 's. 

proof: 

Using lemma 16.11 it is easy to see that Xa^„ fulfills the hypothesis of lemma IA.3I in 
the appendix |Xl and this proves the theorem. □ 

I.e. we have proved the conjecture of Brezin and Hikami [3J. In fact, we notice that 
the property of being a polynomial in the r's, is not specific to angular integrals, but 
comes only from the global symmetries. 

6.1 Recursion without residues for /3 integer 

For /3 G N, the residues in recursion eq. (16- ip can be performed, and they compute 
derivatives of the integrand. Thus eq. (16- ip can be rewritten: 

= AiY^rllix, - x^) Res . . . ^^^^ 

fj^ A,-.a, (Ai - ai)P (A„_i - a„_i)'^ 



/3 



A(r„)2/5 fdV-'^ dx, dx, 



Res 



riLi nr=i,^fc(i/fc - 



rn TTn—l 

1. 

(6-7) 

i.e. we can perform the residues: 



ai=yi 



(6- 
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More explicitly 



TT y^i,nyn,i) "TT 

1/1 d 



(6-9) 



6.2 n = 3 

For 77, = 3, and arbitrary f3, we prove that: 
Theorem 6.2 




(6-10) 

In fact, for (3 integer, the sum over k is finite and reduces to k < (3 — 1. 

The proof, rather technical, is given in appendix O We used the Calogero equation. 



6.3 Conjecture for higher n 

Applying the recursion relations of this article, we also computed the n = 4 case for 
small values of jS: 



I. 



2,4 



and: 



n'52,o(^ii) 



i<j 



(6-11) 



+ TF <52,o(n,2)Q2,o(n,3)Q2,o(n,4) <52,i(t2,3)(52,i(7-2,4)<52,i(t3,4) + Sym 

16 

+ Y^<52,o(n,2) Q2,l(n,3)<52,l(n,4)Q2,l(T-2,3)Q2,l(7-2,4)Q2,l(T-3,4) + Sym 



^3,4 = 64j]Q3,o(r.,) 



i<j 
4 

+ g Q3,o(n,2)<53,o(n,3)<53,o(T-l,4) <53,1 (t-2,3)<53,1 (t-2,4)<53,1 (7-3,4) + Sym 

+ Q3,o(n,2)Q3,o(n,3)Q3,o(n,4) <53,2(t2,3)Q3,2(t2,4)Q3,2(t3,4) + Sym 

+ ... (6-12) 

Those expressions lead us to conjecture a general form in terms of Bessel polynomials: 
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Conjecture 6.1 We conjecture that for all n and (3, Ip^n is of the form: 



{1} 



(6-13) 

Unfortunately, we have not been able so far to determine the general form of the 
coefficients A^iy for n > 3 (except n = 4 and [3 = 2). 

6.4 Symplectic case (5 = 2 

For P = 2, the recursion eq. (l6-8p . reduces to: 



n— 1 



2,n 



i=l 



nLinr=iU(^^- 



n-1 

n 



(6 - 14) 



1=1 k=l,=f=i 



2 d 



ai=yi 



From a recursion hypothesis, we assume that X2,„_i(X„_i, a) is a polynomial in the r's 
of the form: 

{1} i<j 

where for every pair we have G {0, 1}. We recall that: 



(6-15) 



\0>r= Q2,oir)=T^ + T 



Thus we may write: 



d 
dtik 



|1 >.= Qi,o(r) = 2r (6-16) 



ti,k = -^{xi - Xk){ai - Qk) (6-17) 



and we define the operators Ci^k acting on functions of the variable Ti^k such that: 



1 1 d 



Cih — — 



Tifc 2 dtik 



(6-18) 



and all derivatives must be eventually computed at ti^k = Ti^k- 

Since our operators act on expressions of the form eq. (16-151) . we need to compute: 



Cjj|0 >— Cij.Q2fi{ji j) — Cij.{T^j + Tij) — - 



(6-19) 
(6-20) 
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And we may also have terms of the form Cij.Cj^i, for which we have: 
Finally we have: 

«— 1 n— 1 ^ n— 1 

^^,n = n n (i + 7- + E ^ ^^.^) ^v^- L (6-22) 

It is more convenient to rewrite this in terms of a Hilbert space with basis |0 >= (^2,0 
and |1 >= (^2,1, and thus: 



n— 1 ^ n— 1 



(6 - 23) 



2,n-l 



where 



We have: 



A^,k = '^n^kCi^k (6-24) 



A|0>=^|1> , A|1>=|1> (6-25) 



and: 



= 2(A,fe - 1) (6-26) 
Unfortunately we have not been able to go further with this formulation. 

6.4.1 Triangle conjecture 

We have seen that there is an operator formalism for computing angular integrals, with 
operators Ajj associated to "edges" (?,j). However, it can be seen for n = 3,4, that 
operator edges appear only in certain combinations, which involve triangles {i,j,k). 
We thus introduce the triangles operators: 

Ti,j,k = '^i,j'^j,k'^k,i (6-27) 
where tTjj = iTj^i is the projector on state ||1 >i,j= Tij: 

(6-28) 

With this notation we have: 

^2,3 = (1 + ^^1,2,3). n |0>^'^- (6-29) 
l<i<i<3 
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where we recall that |0 >i,j= t^j + Xjj. 
And 

^2,4 = ~'~ 2 ^"^^'^'^ ~'~ "^l'^'*^ ~'~ -^1=3,4 + ^2,3,4) 

+ ^(^1,2,3^1,2,4 + T^i, 2,3^1,3,4 + 7'i,2,47'l,3,4 

+Tl,2,3T2,3,4 + Ti,2,4T2,3,4 + Ti,3,4T2,3,4)). JJ 1^ >M (6-30) 

l<j<jr<4 

We are naturally led to conjecture that: 

^2,n = ( Yi n ) • n io >^.i (6-31) 

triangulations T TeT i<j 

We have: 

= I , C(ij^k) = - , C{i^j^k),{i,j,i) = ^ 5 • • • (6-32) 

We have not been able so far to prove this conjecture. It is to be noted from the 
low values of n, that triangles seem to play a role for all /3 G N. 

6.4.2 Additional results: determinant al recursion 

Just for completeness, we give another form of the recursion eq. (16-141) . in terms of 
determinants: 

det X - Xn+i - y^;^^ + B + dy 

-'2,n+ll^n+li l) — \ l2,n{Xn'i Yn) (6-33) 

det (A - Xn+l) 

where 

/2,„(X, Y) = ^ , J2,„(X, Y) (6-34) 
and where B is the antisymmetric matrix 

B,j = (6-35) 

Vi - Vj 

and dy = diag(9j;^, . . . , dy„). This is proved by observing that the expansion of the de- 
terminant eq. fl6-33p can be interpreted like a Wick's expansion equivalent to eq. (16-141) . 

7 Conclusion 

In this article, we have found many new relations and new representations of angular 
integrals. 
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First, we have been able to rewrite angular integrals with a complicated Haar 
measure, in terms of usual Lebesgue measure contour integrals. Then, we have deduced 
duality and recursion formulae. 

This allowed us to prove Brezin-Hikami's conjecture, and to find some explicit form 
for n = 3, and conjecture some explicit form in terms of Bessel polynomials for the 
general case. 

For f3 = 2, we have simplified our recursion (computed the residues). The same 
method seems to be applicable for higher (3, but we have not done it in this article. 

We have obtained many new forms of angular integrals, but unfortunately, this does 
not seem to be the end of the story. Our expressions are still not explicit enough to be 
useful for computing matrix integrals. The form of our expressions, strongly suggest 
that the kernel determinantal formulae [25] in the /? = 1 case, could be replaced by 
hyperdeterminantal formulae for higher /3, but this is still to be understood. The best 
thing, would be to get expressions with enough structure to generalize the method over 
integration of matrix variables of Mehta |26j . 
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A Appendix: Polynomials of r 

Lemma A.l Let 

Pn{X, Y) = X1X2 . . . Xn-lXn ?/„+l?/n.+2 • • • ?/2n-l?/2n+l/ll/2 • • • yn-lVn Xn+lXn+2 ■ ■ ■ X2n-\X2n 

(1-1) 

We prove that P„, is a polynomial of degree n in the r's, where Tj.2n+i = —\xiyi and 
= —j{xi — Xj){yi — yj), with integer coefficients: 

Pn e Z[r] (1-2) 

proof: 
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It clearly holds for n = and n — 1. Indeed the n — 1 case reads: 

xm + X2yi = xini + X2y2 - {xi - X2){yi - 2/2) = 2(ti,2 - ri,3 - T2fi) (1-3) 

Assume that the lemma holds up to n — 1, and let us prove it for n. In the following 
A = B means that A — B is & polynomial in the r's. 

Pn 



X1X2 . . . Xn-lXn yn+iyn+2 ■ ■ ■ y2n-iy2n 
+I/1I/2 • • • yn-iy-n Xn+lXn+2 ■ ■ ■ X2n-\X2n 



{XiUn- 


K + Xn+iyi){x2. . .Xn 


— iXn 


yn+2- 


+y2 ■ ■ 


■ yn-iy-n Xn+2 ■ ■ ■ X2n- 


-lX2n 


) 


-X2 . 


■ Xn-lXnXn+1 yn+2 ■ ■ 


■y2n- 


-l2/2n2/l 


-y2-- 


■ yn-iynVn+l Xn+2 ■ ■ ■ 


X2n- 


lX2nXl 


-X2 . 


■ Xn—lXnXn+1 yn+2 ■ ■ 


■y2n- 


-mnyi 


-2/2 •• 


■ yn-iynyn+l Xn+2 ■ ■ ■ 


X2n- 


lX2nXl 



(1-4) 

and then: 

Pn 

= —X2 ■ ■ ■ Xn-lXnXn+l yn+2 ■ ■ ■ y2n-iy2nyi 

-y2 ■ ■ ■ yn-iynVn+l Xn+2 ■ ■ ■ X2n-\X2nX\ 
= -{X2yi + Xiy2) (2:3 •• • Xn-lXnXn+1 yn+2 ■ ■ ■ 2/2n-l2/2n 

+2/3 • • • yn-Wnyn+l Xn+2 ■ ■ ■ X2n-\X2n) 

+Xn+\XxXz . . . Xn-lXn 2/22/n+2 • • • 2/2n + yn+Wm ■ ■ ■ yn X2Xn+2 ■ ■ ■ X2n 

(1-5) 

Repeating the same operation recursively we obtain \/k: 

Pn = Xn+l ■ ■ ■ Xn+k Xl Xk+2 • • • a^n 2/2 • • • 2/A:+l 2/fc+n+l • • • 2/2n 

+2/n+l • • • yn+k 2/1 yk+2 ■ ■ - yn X2 ■ ■ ■ Xk+l Xk+n+l ■ ■ ■ X2n 

(1-6) 

In particular for A; = n — 2 we find: 

Pn = Xn+l . . .X2n-2XlXn y2 ■ ■ ■yn-iy2n-iy2n 

+2/n+l • • • 2/2n-2 2/l 2/n 3:2 • • • Xn-l X2n-\X2n 
= {Xiy2n-1 + X2n-iyi)iXn+l ■ ■ ■ X2n-2 2/2 • • • 2/n-l 2/2n 

+2/„+l . . . 2/2n-2 2/n X2 . . . Xn-l X2n) 

-Xn . . . X2n-1 2/l2/2 • • • 2/n-l 2/2n " 2/n • • • 2/2n-l X1X2 • • • Xn-l X2n) 
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(1-7) 



and we repeat the same operations: 

Pn = -Xn ■ ■ ■ X2n-1 2/12/2 ■ ■ ■ Vn-l y2n - Vn ■ ■ ■ y2n-l X1X2 ■ ■ ■ Xn-1 X2n) 

= - {XnUl + XiUn) {Xn+l ■ ■ ■ X^n-l 2/2 ■ ■ ■ Z/n-l y2n 

+yn+l ■ ■ ■ y2n-l X2--- X^-l X2n) 

+Xi X„+i . . . X2n-l 2/2 ■ ■ ■ Z/n y2n + VlVn+l ■ ■ ■ y2n-l X2 ■ ■ ■ Xn X2n 

(1-8) 

and once more 

Pn = XiXn+l---X2n-iy2---yny2n + yiyn+l---y2n-l X2...XnX2n 

= {Xiy2n + X2„yi) {Xn+1 ■ ■ ■ X2n~l 2/2 ■ ■ ■ Z/n + Z/n+l ■ ■ ■ y2n-l X2 ■ ■ ■ Xn) 

-Xn+l ...X2n 2/12/2 ■■■yn- Z/n+l ■ ■ ■ y2n Xi . . . Xn 

(1-9) 

Therefore P„ = —Pm i-e. P„ is a polynomial in the r's. 
□ 

Lemma A. 2 Let 

n n 

PaAX: Y)^ll Xf^yf + n ^^i'yr , J] = J] A = (1-10) 

i=l i=l i i 

We prove that P^^p is a polynomial of degree d in the r 's, with integer coefficients: 

Pn e Z[t] (1-11) 

proof: 

We proceed by recursion on the total degree d = = Pi. The lemma clearly 

holds for d = and d = 1. 

Assume d>2 and that the lemma holds up to c? — 1. We will prove it for d. 

• if there exists i such that ccj/^j > 0, then P^^p is the factor of Ti^n+i times P^',/?' of 
smaller degree, and from the recursion hypothesis it holds. 

• assume that \/i,ai(3i — 0. Since d>2, there must exist some i such that Oii > 1 
and some j such that Pj > I- Let us choose k and I such that: 

ak = max{aj} > 1 , Pi = max{/3j} > 1 (1-12) 

i i 

We have k I. 
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If Ofc = /3; = 1, then we can apply lemma lX!Tl and thus the lemma is proved for 
that case. 

• therefore, we now assume that akPi > 2, and with no loss of generality we may 
assume that > 2. We write: 

Pn 

— ■ ■ - -^n yi ■ ■ ■ y-n ~T~ -^l ■ ■ ■•^n iJl ■ ■ ■ y-n 

{xuyi + xm){xT ■ . . xt'-' ...xl-yt\.. yt^ . . . 2/^ 

"ryi ■ ■ ■ i/fc • • • i/n "''I •••J-; ■■■Xn) 

. . . Xj!. X; . . . yi . . . yi y^ ■ ■ ■ iJn 
yi ■ ■ ■ (Jk yi ■■■yn xi ...Xi •••x„ 
= -Xkyk [x^^ ...x^'^ Xi' ... yl... yp Uk ■ ■ ■ yn 



"Ti/l • • • i/fc ill ■■■yn Xi ...Xi X^i. ...X.„ 



;i - 13) 



From the recursion hypothesis, the RHS is a polynomial in the r's, and thus we have 
proved the lemma. 

□ 

Lemma A. 3 Let P he a polynomial of2n+2 variables xi, . . . , x„, x„+i, yi, . . . , yn, yn+i, 
with the following properties: 

• P is invariant by translations Vi, Xj — Xj + Sx, yi ^ yi + Sy, 

• P is invariant under \/i, Xj — Axj, yi jyi , 

• P is symmetric in the exchange X *-^Y , 
Then P is a polynomial of the Tij 's. 

proof: 

Because of invariance by translation, we can always assume that x„+i = yn+i = 0. 
Then, the other properties imply that P is a linear combination of monomials of the 
type PaA^^ ^) of LemmalAJl 

□ 



B Appendix: Loop equations 

Loop equations for matrix models have been studied for a long time ^7\. Loop equa- 
tions (sometimes called Ward identities or Schwinger-Dyson equations) for /3 ensembles 
can be found for instance in [121 ESI [HI [20l [211 [£]• Here, we summarize the method. 
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• 1-matrix model in eigenvalue representation, for arbitrary 
Consider the integral: 



Z = JdX,...dX^ e- HiXj - Xiff" (2-1) 



i<j 

If we make an infinitesimal local change of variable Aj — > Aj + eX^ + O(e^), we find: 
Z = {l + 0{e')) fdX,...dXnllil + ekXt')e-^i''^'^^ll{l-eX^V'iXi)) 

(2-2) 



i.e., by considering the term linear in e: 

ik-l 



= ( E + /5 E E - E ^'^'(^^) > (2-3) 

\ i 1=0 i^j i I 

It can be written collectively by summing over k with this is equivalent to 

consider a local change of variable Aj — > AjH — ^-l-O(e^), and we write uix) — Y\- 

= (|-u;'(x) + /3(u;(x)^ + u\x)) - ^ (2-4) 



I.e. 



\ x-M/ ^\ x-M x-M/ ^^ ^ dx \ x-M/ ^ ' 



• 2-matrix model for /3 = 1/2, 1, 2: 
Similarly if we consider a 2-matrix model: 



j dM^dM2 e-'^mM,)+v,{M,)-M,M,) (2-6) 

again we make a local change of variable Mi Mi + e^r^g^A + 0(e^). The Jacobian of 
this change of variable is computed as a split rule (cf []), it can be computed for each 
of the 3 ensembles {3 = 1/2, 1, 2 and is worth: 

1^,1 <9 _ , 1 



dMi^dMi ( l + e/?Tr ^TrA 1_ + e(/3 - 1) IVA t— + 0{e^] 

\ x-Mi X- Ml ^ ' dx X- Ml ^ ' 

(2-7) 

Thus we find: 

= -(tv(V/(M0-M.)^)+,5(t.-^T.^ ' 



x-Mi x-Mi 
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+ (ttA^—-) (2-8) 



dx \ X — Ml 

• Without a definition of a 2-matrix integral for arbitrary /?, it is not possible to 
find the loop equation for any [3. However, we see that equation eq. fl2-8p is valid for 
the 2-matrix model for (3 = 1/2,1,2, and is valid for the 1-matrix eigenvalue model 
for any j3. Therefore it is natural to take it as a definition of the 2-matrix model for 
arbitrary (3. 

C Appendix: Proof n = 3 
C.l BESSEL ZOOLOGY 

We define the functions 

«.w-EH'^^^- (3-1) 

We have 



2x'Q',{x)-2PxQ,{x) = 2 



r(/3-0 (/-I) 

-I 



x' Q; (x) -2Px Q', (x) + 2/3 Q, (x) = - ^ ^l^^l^l] |- x^-^ (3-2) 

Changing / into / -|- 1 in 13-21 we obtain the differential equation 

x^ Ql (x) -2x{(3- x) Q'p (x) +2/5(1- x) Qp (x) = (3-3) 

We define the functions 

Q,A.) = (-2)'.'^-' (x4V (3-4) 



^ dx J \ x^ / 

Qpfiix) = Qp{x) (3-5) 

that is 

From [331 we obtain the recurrence 

Qp^k {x) = 2{(3-k + l) Qp^k-i (x) - 2x Q'p^,_, (x) (3-7) 

For instance 

Q^,i (x) = 2f3 Q^,o (x) - 2x Q'^^^ (x) (3-8) 
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and 

Q;3,2 (x) = 2iP-l) Q^,i (x) - 2x Q'f,^, (x) (3-9) 

Qfs,2 (x) = Ax^ Ql, (x) -8xiP- 1) Q'f,^, (x) + 4/5 (/3 - 1) Q^,o (x) (3-10) 

Thus, we have, from 13-81 and 13-10^ the expressions Q'^ q (x) and Q^^o (^) terms of 
(5/3,1 (x) and (5/3,2 (a^)- Equation 13-31 becomes 

Q/3,2 (x) + 4 (1 - x) g^,i (x) - 4/5 (/? - 1) Q/3,0 (x) = (3-11) 

From 13-11^ by derivatives ^ and recurrence, it is easy to show that 

(x) + 4 (A; + 1 - x) Qfs,k+i (x) - 4 [/3 (/? - 1) - A; (A; + 1)] Qfs,k (x) = (3-12) 



C.2 CALOGERO N=3 

We consider the Calogero differential operator 
and we look for solutions 

Hc^lagero ^{Xi, Vi) = (^y^^ ^(Xi, Z/i) (3-14) 

where the solutions $ (xj, ?/i)have a certain number of symmetry properties described 
somewhere else. We write 

<^{X^, y^) = f{x^, Vi) e^^^^^V^ (3-15) 

SO that the equations 13-131 and 13-141 become 

^ ^ ( ^ + ' • ^ "'"^•^ 

doc ^ X 7 ^ V dx i doc ij J dx 7 

i=l « i<j ' ^ * ^ j=l * 

and 

^/(a;., y^)=0 (3-17) 
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Let us introduce the variables 



1 


X2) Y12 


(3-18) 


b = 77 (a^i — a^s) iui — Us) = (xi — 




(3-19) 


c = ^ (3^2 - 2:3) (^2 - 2/3) = (a;2 - 


X3) Y23 


(3-20) 


where 

Y — — y 

and we look for solutions of the type / (a, b, c) . We have 




(3-21) 


and 




(3-22) 
(3-23) 
(3-24) 



dxl 
dxl 



Y^ifa-i + 2F21^23/ac + ^23/c2 
Y^lfy^ + 2Y3iy32/f,c + ^32 /c2 



The equations 13-161 and 13-171 become 

D f (a, b, c) = Di f (a, b, c) + D2 f (a, b, c) 







D,f{a,b,c) = Fj 



12 



2/:.+4(^+i)/:+4^/ 

a a 



D2fia,b,c) = 2Y^2Yi 



13 



/ab + ^ /b + ^ /a 



+ ctrc.perm. 
+ circ.perm. 



We now try the functions 

/ (a, c) = 



^fc(a) Q(3,k{b) Qf3,k{c) {k,k,k} 



a2/3 62/3 ^^^^^2/3 

where ^ (x) are defined in (13-4113-51) and l3-6[ We consider 

Qf3,k (a) 



/(«) 



,2/3 



3-25) 
3-26) 
3-27) 

3-28) 
3-29) 

3-30) 



3-31) 



3-32) 
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we have 



9 f'Mod 2a Q'p,k («)-2/^ Qp,k (g) 



2/3+1 



By derivation we obtain 

2/a2 +'^— fa~ 2^ / = 2,3+2 

so that 

a 



Finally we obtain 

2/:2+4(^ + l)/: + 4^/ = -^ 
a a 2a 



2k Qp^k {a) + Q/3,fc+i (a) 

^2/3+1 



iQ^,fc+2 (a) + (/3 + 2A; + 2) Qp^k+i (a) 
+2A; (/? + A; + 1) g^,fc (a) 



\Ql3,k+2 (a) + 2 (A; + 1) Q^fc+i (a) 
+2A; (A; + 1) Q^,^ (a) 



,2,3+2 



2,3+2 



^Q/3,fc+2 (a) + 2 (/c + 1) (a) 
-2[/?(/3-l)-A;(A;+l)]Q^,fc (a) 



2,3+2 



Q/3,fc+2 (a) + 4 (/c + 1 - a) Q/3,fc+i (a) 
-4[/3(/3-l)-A;(A; + l-2a)]Q^,fc (a) 



Now, we use the recurrence relation 13-121 and get the simple result 

4A; 



2/:2+4(^ + l)/: + 4^/ 



,2,3+1 



■,k (a) 



We just proved that 



(abc) 



2/3 



I b c J [ahcf 



We further transform the result l3-39[ We have 

Yi2 _ Yi2 

{xi - X2) A (a;) 



y2 
-'12 



(Xi - 0:3) {X2 - X3) 



y2 
-'12 



12 



[F - (a:i - 0:2)'] 



(3-33) 
(3-34) 

3-35) 
(3-36) 

(3-37) 

(3-38) 

(3-39) 

(3-40) 
(3-41) 



where 



A (x) = (Xi - X2) {Xi - X3) (X2 - X3) 

F = -|- ^2 + Xg — X1X2 — X1X3 — X2X3 



By circular permutation we also have 



y2 



y2 
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[F - (X3 - x,y] 

A fx) - - ^^^T 



A (x 

^23 



(3-42) 
(3-43) 



(3-44) 
(3-45) 
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We note that in 13-391 the quantity F disappear since 



Yl2 + 1^23 + isi = 



(3-46) 



We may write now 



Ak . \k,k,k} 

[[xi — X2) a + circ.perm.\ 



(abc) 



2/3 



(3-47) 



We now consider 



Using 



we have 



{abcf" 
f{a,b,c) = f{a) f{h) f{c) 



/(b) 



but in 13-301 the term ^ /(a) /(6) /(c) disappears since 



ah 



ah 

Then, from 13-33) we get 



+ circ.perm. = — h circ.perm. = U 



A(x) 



V Q/3,fc(a)+|Q/3,fc+i(a)] [fc Qf(,fc(6)+5Q/3,fc+i(fe)] Qf, fc(c) 
^-'12-'13 



+drc.perm. 



7W 



(3-48) 
(3-49) 

/(c) 
(3-50) 

(3-51) 



(3-52) 



Again, the term containing Q^^^ (o) Qp^k ip) Qp^k (c) disappears by 13-511 We write 



1 



[ahcf^ A (x) (a6c)"^ 

+circ.perm 



{X2 - X3) 



i{fc + l,A; + l,A;} 
+k {k + l,k,k} + k {fc, fc + l,fc} 

(3-53) 



We note that 



(x2 — Xs) [{k + 1, fc, fc} + {k^k + 1, A;} + {fc, k,k + 1}] + circ.perm. = 
so that 



1^2 



1 1 



(abcf^ A (a;) (abc)^" _ 



{X2 - X3) 



{k + l,k + l,k} 

-k {k, k,k + l} 



+ circ.perm. 



(3-54) 



(3-55) 
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We now collect Di and D2. From [3^ 13^ 13^ and 13^ we obtain 



{k,k,k] _ 1 1 
{ahcf^ ~ A (x) (ahcf^ 



D 

Again, 
so that 

j^ {k,k,k} _ 1 1__ 

We now use equation (10) and write 
^ {k,k,k] _ 1 1 



{x2 - a^s) 



Akc {k, k,k} + \{k + l,k + l, k} 
-k {k,k,k + l} 
+circ.perm. 



{x2 — X3) {k, k, k} + circ.perm. = 



(3-56) 
(3-57) 



{X2 - X3) 



Ak {c-k){k,k,k} 
+ ^{k + l,k + l,k} - k {k, k,k + l} 
+circ.perm. 



(abcf" A (x) (abcf" 

Consequently, we obtain the remarquable result 
T{l3-k) 1 {k,k,k} 



(X2 - X3) [i {fc + 1, + 1, fc} 
Ak i/3-k) {l3 + k-l){k,k,k-l}] 
+circ.perm. 



(3-58) 



(3-59) 



D 



T{(3 + k)8''k\ (abcf^ 



1 1 



2A (x) (a6c)'^ 
Clearly enough, we define for (3 not integer 



(X2 - X3) 



r(/3+fc-l)8'=-i (fc-l)! l'^'"'''^ ^/ 

-(-czrc.perm. 



(3-60) 



then, 



r (/5 + A;) 8'= A;! (afec)^^ 
D / (a6c) = + circ.perm. 



Now, if /? is an integer 



Ql3,k>p (x) = 



and we define 



so that 



^r(/? + fc)8'=fc! (afec)^'^ 
D f (abc) = + circ.perm. 



(3-61) 

(3-62) 
(3-63) 

(3-64) 
(3-65) 
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We proved that a solution to 13-141 is 



'oo or 0—1 



-f^Calogero^ {Xi, Vi) = 

where -ffcaiogero is given in 13-131 



T{P-k) 1 {k,k,k} 



^ T{P + k)8'^k\ (abc) 



2/3 



gELi (3.66) 
(3-67) 
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